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The Special Affine Fourier Transform or the SAFT generalizes a number of well known unitary 
transformations as well as signal processing and optics related mathematical operations. Unlike the 
Fourier transform, the SAFT does not work well with the standard convolution operation. 

Recently, Q. Xiang and K. Y. Qin introduced a new convolution operation that is more suitable 
for the SAFT and by which the SAFT of the convolution of two functions is the product of their 
SAFTs and a phase factor. However, their convolution structure does not work well with the inverse 
transform in sofar as the inverse transform of the product of two functions is not equal to the 
convolution of the transforms. In this article we introduce a new convolution operation that works 
well with both the SAFT and its inverse leading to an analogue of the convolution and product 
formulas for the Fourier transform. Furthermore, we introduce a second convolution operation that 
leads to the elimination of the phase factor in the convolution formula obtained by Q. Xiang and 
K. Y. Qin. 
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INTRODUCTION 


Let £1' :/—>■/ be a unitary, integral operator which maps / to its transform domain representation /. For example, 
if is the Fourier operator, that is, 

^FT [/] (^) = = ^ [ f it) (1) 

VJtf Jr 

then / is identified as the frequency domain of /. Furthermore, let * denote the standard convolution operator defined 

by 

if*9){t) = -^ [ f{x)g{t-x)dx. 

VJtt Jr 

It is well known that for some functions / and g, and 17 = JJft, we have, 

•^FT [/ * g] (w) = [/] (w) 31ft [ 5 ] (‘^) ■ (2) 

This result is known as the Fourier convolution theorem. 

The purpose of this paper is to extend and establish the Fourier convolution theorem for the Special Affine 
Fourier Transform (SAFT)—a phase space transform—which generalizes a number of well known transformations. 



2 


Technical Report No. XX, Vol. XX, No. X, Month 20XX 


KAs (t, u) 


exp + 2t{j) — uj) — 2uj [dp — bq))^ , 


Kb = 



(3) 


Some of the interesting transformations and signal/optical operations that can be obtained from the SAFT as special 
cases are listed in Table I. 


Phase Space Transformations 

Phase-Space transformations such the the fractional Fourier Transform (FrFT) [1] and the Linear Canonical Trans¬ 
form (LCT) [2] are becoming increasing popular in the areas of signal processing and communications. A remarkable 
feature of the phase space transformations is that they generalize the Fourier Transformation and hence, all the 
mathematical developments are compatible with the Fourier analysis. 

In recent years, a number of fundamental, signal processing centric theories for phase space have been developed. 
Some examples include convolution theorems [3-8], sampling theory [8-14], time-frequency representations [15, 16], 
shift-invariant signal approximation [8], sparse sampling theory [17] and super-resolution theory [18]. 


Convolution Theorems for the PrFT and the LCT 

In context of signal processing theory, Almeida first studied the fractional Fourier Transform (FrFT) domain 
representation of convolution and product operators [3]. Unfortunately, Almeida’s formulation did not conform with 
the classical Fourier convolution-multiplication property. That is to say, the convolution of functions in time domain 
did not result in multiplication of their respective FrFT spectrums. As a follow up, Zayed formulated the convolution 
operation for the FrFT which resulted in an elegant convolution-multiplication property in FrFT domain [4]. 

Recently, Xiang and Qin [7] introduced a new convolution operation that is more suitable for the SAFT and by 
which the SAFT of the convolution of two functions is the product of their SAFTs and a phase factor. However, their 
convolution structure does not work well with the inverse transform in sofar as the inverse transform of the product 
of two functions is not equal to the convolution of the transforms. 

In this article we introduce a new convolution operation that works well with both the SAFT and its inverse leading 
to an analogue of the convolution and product formulas for the Fourier transform. Furthermore, we introduce a second 
convolution operation that leads to the elimination of the phase factor in the convolution formula obtained in [7[. 


Special Affine Fourier Transform (SAFT) 

The SAFT was introduced by Abe and Sheridan [19, 20] who studied a transformation in phase space that was 
associated with a general, inhomogeneous, lossless linear mapping. Such a transformation can model a number of 
optical operations such as rotation and magnification (see Table I). 

Let f* denote the complex-conjugate of / and (/, g) = J f (t) g* (f) dt be 
SAFT operation, that is, '■ f —> /as, is defined as, 

? ( \-\ (•>‘^)) ^ 

(d(a; — p)) b 

where, 

f 2 X 3 

• Ag ' is the augmented SAFT parameter matrix of form. 

As = [ A I A ] (5) 

which is in turn parameterized by the LCT matrix Al [2[ (see Table I) and an offset vector A such that, 

A = [ “ ^ ] with ad — be = 1 and A = [ g ] ■ 

This is the reason the SAFT is sometimes referred to as the Ojfset Linear Canonical Transform or the OLCT. 


the standard L 2 inner-product. The 


^0 
= 0 ’ 


(4) 
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• i^As in (3) is the SAFT kernel parameterized by SAFT matrix As- 

Thanks to the additive property of the SAFT/OLCT [21], the inverse-SAFT (or the iSAFT) is simply the SAFT 
evaluated using matrix A™'' with parameters, 


'+d -b 

bq- 

dp 


'+d -b 

Po 

—c +a 

cp — 

aq 


—c +a 

qo 


As a result, we are able to define the inverse transform iSAFT, 

/ (0 = C'a'- (/as. (•-O) 


( 6 ) 

(7) 


with some transform dependent phase constant. 


CAmv = exp [cdp^ + ahq^ — 2adpq)^ . 


Next, we develop a convolution structure for the SAFT denoted by *As so that we can obtain a representation of 
form. 


•^AFT [/ *As g] oc =3^AFT [/] -^^AFT [ 5 ] 
which is consistent with the Fourier convolution theorem (2). 


CONVOLUTION THEOREM FOR THE SAFT 


Before we define the convolution operation in the SAFT domain, let us introduce the chirp modulation operation. 
Definition 1 (Chirp Modulation). Let A = [aj^k\ he a 2 x 2 matrix. We define the modulation function, 

rriA {t) exp • (8) 

Furthermore, for a given function /, we define its chirp modulated functions associated with the matrix A as, 

7 {t) ruA (t) f (t) and J (t) ‘^= m\ (t) f (t). (9) 


For example, let A = As, then, we have, f (t) = mj\_^f {t) = f {t). For the case when A = A™'', we get, 

7 (t) = ruAinv {t) f (t) = e-^'^fit). 

Next we define the SAFT convolution operator. 

Definition 2 (SAFT Convolution/Filtering). Let f and g he two given functions and * denote the usual convolution 
operation (see (2)). The SAFT convolution is defined hy, 

h (t) = (/ *As g) (t) ATfemAj (t) (j (t) * g (f)^ . (10) 

Figure 1 illustrates the block diagram for SAFT domain convolution defined in Definition 2. Next, we state the 
convolution and product theorem for the SAFT domain. 

Theorem 1 (SAFT Convolution and Product Theorem). Let f and g he any two given functions for which the 
convolution *As exists and set. 


h {t) = if *As g) (t) ■ 

Furthermore, let /as (w) ,gAs {^) ^As (w) he the SAFT of f,g and h, respectively. Then we have, 

h it) = if *As g) it) hAs (w) = $As (w) 7as (w) 5As (w) , 
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Kte- 


f - 


e-> 2b 


4 


9 -Mg)— 9 


Convolution 



/ *As 9 


FIG. 1. Conceptual definition of SAFT convolution. 


where (w) = . Moreover, let, 

h (t) = $A*™ (t) f (t) g {t) with $Aj™ (t) = 

then, we have, hAs (w) = Ca^™ (j *Ai™ (w). 

Proof. We begin with computing the SAFT of h, 


(4) 


hAs = '^SAFT [h] (w) = {h (t) , KAs (i, w)) 

= I h{t) K* {t, oj) dt 

Js. 

( 10 ) 


Kbm*A (t) / f {z) g {t- z) dz k* {t, oj) dt 




WAs (t) rriAs (t) / f{z)g{t- z) dz dt 


C{i^) 


In the above development, note that the items in the box cancel one another because m\^mAs = 1 (see Definition 1). 
Setting t — z = V and using (9), we obtain an integral of separable form, that is, /ias = If (w) Ig (oj) because, 

Kb^As (^) [[ f i^) "4As (x) g (v) niAs (x) e^^‘^P~‘^'>dxdv 

J Jr 

/lAj (bj) 

= If{‘^)Ig(‘^) = hAs{i^), (11) 

where, for a given function /, we define. 

If (oj) =*' (oj) J f (z) niAs (z) e^i^P~‘^'^dz. (12) 

Indeed, using (12) and (4), it is easy to see that. 

If (oj) = (uj) 7as (w) (13) 

and this result extends to Ig (w) by symmetry. 

We conclude, 

hAs (^) If(c0)lg(cj) 

Y^^As (w)$As (w) 7as (w) • (‘^)^As (w)?As (w) 

= $As (w) /as (w) ?As (w) 


which is the statement of part I of Theorem 1. 

Now we establish the product theorem for the SAFT, 
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Since the inverse-SAFT is the SAFT of a function with As = in (6), we have, 

h{t) 

/ hAs (w)e-^-55- 

Jr 

(i) / hAs{oj)e-^^e-^^^dco. 


By setting, 


/lAs (w) = C'A.nv (J 9 ) (w) 

= C'A.nv (w) (w) * g (w)^ , 

where, for As = A™'' (see (6)), we have, 

_ d()5 ^ 

"5a|"v (•) = e ^ 26 and / (w) = TOAi,nv (w) / (w). 

Upon simplification, we obtain the separable integrals, 

h (t) = (CA^nv Aft*) (t) f f f (^) "lA-- (^) gi^^-v) mAinv (w - 12 ) dve~^ ’ dw 






= y / (z/) e ^ J 9 (^) e ^ ^ dto 

= If (t) 4 (0, 


where. 


4 (t) = CA.nv Aft* it) [ 7 (cc)e-^'^da;. 

V s Jjj 

= W/W- 

As a result, we have, 

h{t)=If{t)Ig{t) 

= (^) / (^) • yn^'J’Ar (^) 9 {t) 

= (i) / (0 9 (t) 

which is the desired result. ■ 


COMPARISON AND ALTERNATIVE RESULTS 

In this section we compare our results with those of [7]. We show that our approach is not only easier to derive, 
but also provides more symmetric formulas to implement. Our convolution formula is the same as the one given in 
[7], both assert that the SAFT of the convolution of two functions is the product of their SAFT and a phase factor 
given by $As (w). But our product formula is different from that in [7] which states that the SAFT of the product of 
two functions / and g is, 

Kb^ls (^) (/As (w) (w) * 9A,r (^) ) . (14) 

The reason our convolution and product formulas are more symmetric and simpler goes back to our definition of 
the chirp modulation, see Definition 1 , which uses the adaptive matrix As that accommodates both the forward and 
backward SAFT. 

Furthermore, we will now derive another convolution for SAFT which eliminates the phase factor <I>As from the 
convolution formula. 
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Definition 3 (Phase-free SAFT Convolution). Let f and g be two given functions and * denote the usual convolution 
operation. The second SAFT convolution * is defined by, 

Ht) = {f*9){t) (t) ^7 (.) * g (.)^ {'/2t). 

In view of this SAFT-convolution, we have the following theorem. 

Theorem 2. Let h{t) = {f ■kg){t). Then, we have, 

hAs (w) = 7 ai (w/?Ai (uj/V2j , 

where /a^ denotes the SAFT of f with respect to the matrix Ai = [A | A/-\/2] (cf. (5)). 

Proof. Let Lt = bq — dp. We have, 

{uj) = V2Kle~'~^ [ 

Jr 

X f f{T)e^^g{V2t — dr. 

Jr 


Setting X = \/2t — r and simplifying the integrals, we obtain, 

hAsiuj) = Kl^X^uj) f e^f{T)dT [ 

Jr Jr 

Jr Jr 

Jr Jr 


But since, 


^A;(2<^)X2C2n^) 


=e 

it follows, that /ias (^) = where 

r ^faP+d(^y+V2pr-V2Tui+V2n{ui/V2)\ 

If[ui) = Kb / e V ^ ^firjdi 

JR 


and similar expression for Ig{uj). But it is easy to see that 

//(w) = 7ai (uj/V2j , 

and this completes the proof. ■ 


Relation to Convolution theory of LCTs: In the special case where p = q = 0 As = Alct, the SAFT reduces 
to the LCT and the last convolution theorem takes the simple form 

■^lct[/ = c31ct[/] (a}/V2j SAcjlg] (w/-\/2^ . 


CONCLUSION 

In this letter, we introduced two definitions of convolution operation that establish the convolution-product theo¬ 
rem for the Special Affine Fourier Transform (SAFT) introduced by Abe and Sheridan [19, 20]. Our result is quite 
general in that our convolution-product theorem is applicable to all the listed unitary transformations in Table I. 
Furthermore, we also presented a product theorem for the SAFT which establishes the fact that the product of 
functions in time amounts to convolution in SAFT domain. We conclude that our construction of the convolution 
structure for the SAFT domain establishes the SAFT duality principle, that is, convolution in one domain amounts 
to multiplication in the transform domain and vice-versa. Our results can be used to develop the semi-discrete 
convolution structure [8] for sampling and approximation theory linked with the Special Affine Fourier Transform, 
but this will be done in a separate project. 
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TABLE I. SAFT, 

Unitary Transformations and Operations 

SAFT Parameters (As) 

Corresponding Unitary Transform 

1 a 6 1 0 1 _ A 

L c d 1 0 J 

Linear Canonical Transform 

r cos 9 sin 9 \ pi _ aO 

L — sin 9 cos 0 1 g J ^ 

Offset Fractional Fourier Transform 

f cos 9 sin 0 101 A 

L - sin 0 cos 0 1 0 1 ^ ■‘*■0 

Fractional Fourier Transform 

r 0 11 pi _ ao 

L - 1 0 1 q\ — ■''■FT 

Offset Fourier Transform (FT) 

[-10 o] = •'^FT 

Fourier Transform (FT) 

[?^!S]=Aft 

Laplace Transform (LT) 

r J cos 0 J sin 0 01 

Ljsin0 — JCOS0 oj 

Fractional Laplace Transform 

r 1 b 1 01 

Lo 11 oJ 

Fresnel Transform 

]b\ 0 1 

Li 1 1 ol 

Bilateral Laplace Transform 

ri -jb 1 01 5 > 0 

Lo 1 1 0 J ’ — ^ 

Gauss—Weierstrass Transform 

1 r 0 1 0 ] 

V2 L 1 0 J 

Bargmann Transform 

SAFT Parameters (As) 

Corresponding Signal Operation 

1 I/O 0 1 0 1 A 

Time Scaling 

ri 0 1 rl _ A 

Lo 1 1 0 J — 

Time Shift 

r 1 0 1 01 _ A 

Lo i| d 

Frequency Shift 

SAFT Parameters (As) 

Corresponding Optical Operation 

n cos 0 sin 0 1 Cn A 

L - sin 0 COS 0 1 0 J ~ 

Rotation 

fl 01 01 _ A 

Lr 1 1 0 J — ■''•T 

Lens Transformation 

ri r? 1 01 _ A 

Lo 11 oJ — 

Free Space Propagation 

L 0 1 0-1 ^ 

Magnification 

r cosh a sinh a 1 01 _ A 

L sinh a. cosh a 1 0 J 

Hyperbolic Transformation 


* ayush@MIT.edu or Ayush.Bhandari@googlemail.com 

[1] H. M. Ozaktas, M. A. Kutay, and Z. Zalevsky, The fractional Fourier transform with applications in optics and signal 
processing. Wiley New York, 2001. 

[2] M. Moshinsky and C. Quesne, “Linear canonical transformations and their unitary representations,” Journal of Mathemat¬ 
ical Physics, vol. 12, no. 8, pp. 1772-1780, 1971. 

[3] L. Almeida, “Product and convolution theorems for the fractional Fourier transform,” IEEE Signal Proc. Lett., vol. 4, no. 1, 
pp. 15-17, Jan 1997. 

[41 A. I. Zayed, “A convolution and product theorem for the fractional Fourier transform,” IEEE Signal Proc. Lett., vol. 5, 
no. 4, pp. 101-103, 1998. 

[5] O. Akay and G. Boudreaux-Bartels, “Fractional convolution and correlation via operator methods and an application to 
detection of linear FM signals,” IEEE Trans, on Sig. Proc., vol. 49, no. 5, pp. 979-993, 2001. 

[6] D. Wei, Q. Ran, Y. Li, J. Ma, and L. Tan, “A convolution and product theorem for the linear canonical transform,” IEEE 
Signal Proc. Lett., vol. 16, no. 10, pp. 853-856, 2009. 

[7] Q. Xiang and K. Qin, “Convolution, correlation, and sampling theorems for the offset linear canonical transform,” Signal, 
Image and Video Processing, pp. 1-10, 2012. 

[8] A. Bhandari and A. 1. Zayed, “Shift-invariant and sampling spaces associated with the fractional Fourier transform domain,” 
IEEE Trans, on Sig. Proc., vol. 60, no. 4, pp. 1627-1637, 2012. 

[9] X.-G. Xia, “On bandlimited signals with fractional Fourier transform,” IEEE Signal Proc. Lett., vol. 3, no. 3, pp. 72-74, 
1996. 

[10] A. Stern, “Sampling of compact signals in offset linear canonical transform domains,” Signal, Image and Video Processing, 
vol. 1, no. 4, pp. 359-367, 2007. 

[11] R. Tao, B.-Z. Li, Y. Wang, and G. K. Aggrey, “On sampling of band-limited signals associated with the linear canonical 
transform,” IEEE Trans, on Sig. Proc., vol. 56, no. 11, pp. 5454-5464, 2008. 

[12] R. Tao, B. Deng, W.-Q. Zhang, and Y. Wang, “Sampling and sampling rate conversion of band limited signals in the 
fractional Fourier transform domain,” IEEE Trans, on Sig. Proc., vol. 56, no. 1, pp. 158-171, 2008. 

[13] H. Zhao, Q.-W. Ran, J. Ma, and L.-Y. Tan, “On bandlimited signals associated with linear canonical transform,” IEEE 
Signal Proc. Lett., vol. 16, no. 5, pp. 343-345, 2009. 

[14] J. Shi, X. Liu, X. Sha, and N. Zhang, “Sampling and reconstruction of signals in function spaces associated with the linear 
canonical transform,” IEEE Trans, on Sig. Proc., vol. 60, no. 11, pp. 6041-6047, 2012. 









Technical Report No. XX, Vol. XX, No. X, Month 20XX 


[15] C. Capus and K. Brown, “Short-time fractional Fourier methods for the time-frequency representation of chirp signals,” 
The Journal of the Acoustical Society of America, vol. 113, no. 6, p. 3253, 2003. 

[16] R. Tao, Y.-L. Li, and Y. Wang, “Short-time fractional Fourier transform and its applications,” IEEE Trans, on Sig. Proc., 
vol. 58, no. 5, pp. 2568-2580, May 2010. 

[17] A. Bhandari and P. Marziliano, “Sampling and reconstruction of sparse signals in fractional Fourier domain,” IEEE Signal 
Proc. Lett., vol. 17, no. 3, pp. 221-224, 2010. 

[18] A. Bhandari, Y. Eldar, and R. Raskar, “Super-resolution in phase space,” in Proc. of IEEE International Conference on 
Acoustics, Speech and Signal Processing (ICASSP), April 2015. 

[19] S. Abe and J. T. Sheridan, “Optical operations on wave functions as the abelian subgroups of the special affine Fourier 
transformation.” Optics Letters, vol. 19, no. 22, p. 1801, Nov 1994. 

[20] S. Abe and J. Sheridan, “Generalization of the fractional Fourier transformation to an arbitrary linear lossless transforma¬ 
tion an operator approach,” J. of Physics A: Math, and General, vol. 27, no. 12, p. 4179, 1994. 

[21] S.-C. Pei and J.-J. Ding, “Eigenfunctions of the offset Eourier, fractional Fourier, and linear canonical transforms.” JOSA 
A, Optics, Image Science, and Vision, vol. 20, no. 3, pp. 522-532, Mar 2003. 


